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We have theoretically analyzed coherent nuclear-spin dynamics induced by electron transport 
through a quantum-dot spin valve. The hyperfine interaction between electron and nuclear spins in a 
quantum dot allows for the transfer of angular momentum from spin-polarized electrons injected from 
ferromagnetic or half-metal leads to the nuclear spin system under a finite voltage bias. Accounting 
for a local nuclear-spin dephasing process prevents the system from becoming stuck in collective 
dark states, allowing a large nuclear polarization to be built up in the long-time limit. After 
reaching a steady state, reversing the voltage bias induces a transient current response as the nuclear 
polarization is reversed. Long-range nuclear-spin coherence leads to a strong enhancement of spin- 
flip transition rates (by an amount proportional to the number of nuclear spins) and is revealed 
by an intense current burst, analogous to superradiant light emission. The crossover to a regime 
with incoherent spin flips occurs on a relatively long time scale, on the order of the single-nuclear- 
spin dephasing time, which can be much longer than the time scale for the superradiant current 
burst. This conclusion is confirmed through a general master equation. For the two limiting regimes 
(coherent/incoherent spin flips) the general master equation recovers our simpler treatment based on 
rate equations, but is also applicable at intermediate dephasing. Throughout this work we assume 
uniform hyperfine couplings, which yield the strongest coherent enhancement. We propose realistic 
strategies, based on isotopic modulation and wavefunction engineering in core-shell nanowires, to 
realize this analytically solvable “box-model” of hyperfine couplings. 

PACS numbers: 68.65.Hb,72.25.Pn,75.78.-n 


I. INTRODUCTION 

Continuous efforts toward the implementation of quan¬ 
tum computing with electron-spin qubits have led to 
several advances in the spin manipulation of single 
electrons in III-V semiconductor quantum dots and 
to improved understanding of the associated decoher¬ 
ence mechanismsP® From these studies, it has become 
clear that hyperfine interactions with nuclear spins in 
the host material typically limit electron spin coher¬ 
ence. Several approaches have been developed to limit 
this source of decoherence (including, e.g., nuclear-spin 
state narrowing through passive measurement 5 ® or ac¬ 
tive feedback controP^Hm 0 f the nuclear Overhauser 
field, or through spin-echo and more general dynamical- 
decoupling techniqueP^ 17 ). These approaches are of¬ 
ten limited to an inconvenient parameter regime (e.g., 
large magnetic field) or require additional pulsed or 
continuous-wave excitation. For these reasons, new 
strategies to accurately control coupled electron-nuclear 
spin dynamics are being actively pursued. 

Quantum dots confining single holes have recently 
emerged as a promising alternative platform for spin 
qubits, since their Ising-like hyperfine coupling allows 
for superior control of coherence, relative to electron 
systemsP^S Another strategy is to exploit group-IV ma¬ 
terials: C, Ge, and Si, all of which can be isotopically en¬ 
riched to be nuclear-spin free! 27 28 However, the nuclear- 
spin bath can also serve as a useful resource, providing a 


highly local tunable effective magnetic field, or long-lived 
quantum memory, as demonstrated by the transfer of 
the qubit state from the electron to the nuclear-spin sys¬ 
tem in NV centers and phosphorus donors in SiPBESJ For 
quantum dots, being able to engineer a well-understood 
form of hyperfine coupling which allows better control 
of the coupled electron-nuclear spin dynamics could lead 
to improved quantum memories that exploit collective 
nuclear-spin degrees of freedonP^l and a fundamentally 
improved level of control over electron-spin coherence. 

A major difficulty in modeling spin dynamics arises 
from the typical non-uniformity of the electron-nuclear 
hyperfine coupling strength. In fact, while many ex¬ 
act and approximate theoretical approaches to strongly- 
coupled electron-nuclear spin dy namics have been de¬ 
veloped for various specific limits)- 5 M33JE2I a controlled 
theory applicable to a typical number of 10 5 -10 6 nu¬ 
clear spins does not exist for many experimentally rel¬ 
evant regimes (e.g., very low magnetic field). On 
the other hand, a simple exact solution based on to¬ 
tal angular-momentum eigenstates can be found if the 
hyperfine co upling-c onstants are uniform, the so-called 
“box-model” ! 34 l 43 l l3 D ue to the simplicity of this limit, 
there is a large relevant body of theoret ical literature, 
addressing electron-spin dephasing) 43 * 45 l 46 J as well as ma¬ 
nipulation and entanglement generation/preservation for 

the electron-spirP^®2 0 r nuclear-spin systemf^^ 

The inhomogeneous hyperfine coupling in current de¬ 
vices is due to the spatial dependence of the electronic 
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wavefunction. For a bound state, this inhomogeneity is 
thus impossible to avoid in III-V materials (where all 
isotopes have a finite nuclear spin). Although uniform 
coupling is often assumed with the purpose of gaining 
insight into realistic setups, it is not always clear to what 
extent these results are applicable to experimentally real¬ 
izable situations. With these motivations in mind, here 
we discuss how Si/Ge core/shell nanowires can offer a 
promising route to realize the uniform limit of the hy- 
perfine interaction. In the following sections, we describe 
quantum-dot designs based on this type of nanowire, 
which approach the uniform-hyperfine-coupling regime in 
a systematic fashion. 

An immediate consequence of uniform “box-model” 
hyperfine couplings would be a strong enhancement of 
the spin-flip electric transport through quantum dots, in¬ 
duced by the flip-flop component of the hyperfine inter¬ 
action. This pheno meno n has been predicted in earlier 
works by Eto et al!— I--J and is analogous to superradi¬ 
ant photon emi ssio n of atoms with uniform coupling to 
the optical field/’ 5 ® An alternative setup with unpolar¬ 
ized contacts has also been recently analyzed in Ref. [57] 
These ideas are also relevant to quantum dots with mag¬ 
netic impurities®® an d, j n addition to transport, opti¬ 
cal superradiance due to the nuclear-spin system has also 
been proposed for single-photon emitters. 1 ^ 

To maximize the influence of spin-flip contributions 
to current, here we will consider a quantum-dot spin 
valve with antiparallel ferromagnetic leads. Experimen¬ 
tal realizations of quantum-dot spin valves include InAs 
quantum dots with Ni ferromagnetic contacts®® as 
well as a spin valve based on a Si nanowirej 65 which 
is compatible with our proposed implementation of uni¬ 
form hyperfine interaction. There is also ongoing inter¬ 
est in carbon nanotubes with ferromagnetic contacts, in¬ 
cluding quantum dot spin valves (see, e.g., the recent 
Ref. l66iand[67]and references therein). With this exper¬ 
imental progress in mind, here we examine a scheme to 
demonstrate superradiance-like behavior in a quantum- 
dot spin valve with ferromagnetic leads. To account for 
imperfectly-polarized ferromagnetic contacts, we allow 
for nuclear-spin states that are both partially polarized 
and fully dephased, as well as spin-flip tunneling pro¬ 
cesses allowing for spin flips in both directions. These 
aspects are not present in the superradiant-like transport 
in an id eal ‘spin -blocked’ regime discussed in previous 
literature 51 52 5 -3 (in our case, for half-metal leads). These 
features are, however, generically relevant for most real¬ 
istic setups. Furthermore, the quantum-dot spin valve 
discussed here offers certain advantages with respect to 
control over the nuclear spin polarization. In particu¬ 
lar, in the setup described here, it is possible to reverse 
the nuclear-spin magnetization direction by reversing the 
bias. This allows for an efficient initialization mechanism 
(without inverting the magnetic field), which is required 
to detect a strong enhancement of the transient current 
to provide evidence of long-range nuclear-spin coherence. 

This paper is organized as follows: in Sec. [ll]we define 


the hyperfine interaction with uniform coupling strength 
and in Sec. Ill we describe possible strategies to approach 
this limit, based on Si/Ge core/shell nanowire quantum 
dots. The rest of the paper analyzes a quantum-dot 
spin valve setup described in Sec. [IVj We introduce rele¬ 
vant tunneling rates in Sec. [V] and discuss the superradi¬ 
ance analogy for half-metal leads in Sec. |Vg In Sec. m\ 
we characterize the spin valve in the two extreme lim¬ 
its of nuclear-spin dephasing (fully coherent/incoherent 
nuclear-spin evolution). We analyze stationary states, 
magnetization dynamics, and the electron current. In 
Sec. |VIIl| the two limiting regimes of Sec. |VII| are obtained 
from a more general master equation, which also allows 
us to discuss the case of intermediate dephasing. Sec¬ 
tion |IX| contains our final remarks and in Appendices [A] 
and |B| we provide some additional details on the exact 
eigenstates with uniform hyperfine coupling and on the 
master equation approach of Sec. 


II. HYPERFINE INTERACTION 


The Fermi contact interaction of a single electron with 
its nuclear bath is well known: 

Hm = A k lk • S, (1) 

k 


where S = | XZ CT a , a^o-id^du^ is the electron spin opera¬ 
tor (we set h = 1), a is the vector of Pauli matrices, d' a 
creates an electron with spin a =7,/. i n the lowest orbital 
of a quantum dot, and Ifc is the spin operator for nucleus 
k. The hyperfine coupling is A k = uo^4 lfc |V ; ( r fc)| 2 j where 
jk indicates the isotopic species, Vq is the atomic volume, 
and ip(rk) is the value of the electronic envelope wave 
function at the nuclear site r k- 

The central spin problem resulting from Eq. 0 is in 
general very complex to analyze!- 4 * 5 * 23 T 33 ^ 8 l 40 * 41 f The limit 
of N nuclei with uniform hyperfine couplings (which al¬ 
lows for an exact solution) is therefore especially inter¬ 
esting: 

H bI = H% + Hy f , (2) 

H “ = | I Z S Z , Hf f = A (. I+ S _ + I-S+), (3) 


where I = Ife is the total nuclear angular momentum, 
I± = I x ± il y , S± = S x ± iS y , and, for later convenience, 
in Eq. ([3]) we have explicitly written the secular and flip- 
flop contribution, H^l and respectively. An exact 

solution for the dynamics of H hf can be found in terms of 
total nuclear angular-momentum eigenstates | I,I Z ) (see, 
e.g., Refs. 03J{5T| and [52] for a solution in the absence of 
a magnetic field and Refs. 03] and 04] for the case with 
a magnetic field), where we omit here and in the follow¬ 
ing an additional permutation quantum number.® The 
underlying assumption of Eq. ([3]) is 


A 3k |V’(rfc)| 2 


-4? for 1 < k < N 

v 0 N — — 

0 otherwise 


( 4 ) 
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FIG. 1. Schematic of a nuclear-spin-free core/shell structure 
with an embedded island of nuclear spins. The core is com¬ 
posed of regions 1-2-3, where the middle section (blue region 
2) is enriched by a spinful isotope. Regions 4 and 5 rep¬ 
resent different shells of the nanowire. By using additional 
gates along the nanowire the electron wavefunction ip( r) is 
centered around region 2, such that the hyperfine couplings 
are approximately uniform. 


Additional electrostatic gates make it possible to lo¬ 
calize a single electron with wavefunction centered at the 
island of nuclear spins, see Fig. [l] Although this does 
not strictly implement Eq. the hyperfine coupling is 
approximately uniform for a wavefunction V’( r ) °f com¬ 
paratively large extent with respect to the size of the 
active central region. A typical core diameter is 10 nm 
and a few atomic layers along the nanowire, with several 
hundred nuclear spins each, can be strongly coupled to a 
bound electron spin. 

In addition to the nearly uniform coupling, a further 
advantage of having an isolated island of nuclear spins is 
that diffusion of nuclear spin polarization out of the quan¬ 
tum dot becomes impossible (in contrast, e.g., to GaAs 
quantum dots). Finally, the electron can be moved away 
from the region 2 through electrostatic gates, turning off 
the contact hyperfine interaction. 


This “box-model” is difficult to implement in actual de¬ 
vices. For the most-studied example of GaAs quantum 
dots, Eq. |4]) becomes impossible to realize, since 7 ^ 0 
for all isotopes of Ga and As and a uniform wavefunc¬ 
tion i/j( r) which abruptly vanishes outside the quantum 
dot is not realistic. On the other hand, it is possible to 
devise strategies approaching Eq. Q in Si/Ge core/shell 
structures, a promising platform for quantum informa¬ 
tion processingPSEI 


B. Wavefunction engineering 

The strategy indicated in the previous section relies 
on confining the nuclear spins to a small region in which 
r) is approximately uniform, but now we describe a 
possible route to realize wavefunctions with large uni¬ 
form regions. In principle, the desired ip{ r ) can always 
be defined through suitable confinement, as can be seen 
for one-dimensional systems. Given a ground-state wave- 
function ip(z) with energy Eq, the corresponding poten¬ 
tial is: 


III. UNIFORM HYPERFINE COUPLING IN 
SI/GE CORE/SHELL NANOWIRES 


V(z) 


1 d 2 i/>(z) 

2 mijj{z) d z 2 


+ A'O; 


( 5 ) 


To realize a uniform hyperfine coupling, it is necessary 
that the host crystal has both isotopes with and without 
nuclear spin, a condition which is indeed satisfied for both 
Ge and Si. We have found that core/shell nanowires offer 
the possibility to approach the ideal situation of uniform 
hyperfine couplings, see Eq. Q. The general structure 
of such a nanowire is illustrated in Fig. |T] with a core 
and multiple shells of different SiGe alloys in the radial 
direction. 72 ^ 


A. Embedded islands of nuclear spins 

A possible design achieving almost uniform hyperfine 
interaction is illustrated in Fig. [l] and is based on in¬ 
cluding nuclear spins only in a small region of the wire 
(region 2). The structure is realis tic to grow, being anal¬ 
ogous to embedded quantum dot spill w here the embed¬ 
ded island is of a different chemical composition than the 
surroundings. Instead, the same alloy can be used here 
for all regions 1-4 (with shell 5 providing confinement), 
but region 2 is isotopically enriched with nuclear spins. 
Isotopic modulation of silicon nanowires has already been 
demonstrated in Ref. m 


with V{z) = Eq constant in the central region where 
ip{z) is uniform. A specific example is simply obtained 
from the following ip{z), defined for — W < z < W and 
uniform in the central region — W + ^ < z < W — 

( sinir{z + W)/d for - W < z <-W + f, 
ip(z) oc< 1 for — W + ^ < z <W — |, 

{ simr(W - z)/d for W - f < z < W. 

( 6 ) 

This wavefunction corresponds to the following potential: 

fOT ~ W < z < —W + f, 

V(z) = i 0 for -W+\ <z<W-$, (7) 

forW-l<z<W, 

where we have have set Eq = 0 in Eq. (|5j). With a differ¬ 
ent choice of ground-state wavefunction, it is also easy to 
obtain from Eq. § a smooth potential V(z) analogous 
to Eq. |7|, i.e., with a central region in which V(z) = 0 
and two negative potential wells on both sides of such a 
central region. 

A similar strategy can be applied to the radial con¬ 
finement and, in fact, core/shell structures allow for the 
realization of a stepwise potential analogous to Eq. ([TJ). 
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FIG. 2. Value of qR ou t — qRin as a function of qRin, as given 
by Eq. (10 1 . The dashed line is the asymptotic value tt/ 2, 
see Eq. m The inset illustrates the radial dependence of 
the confining potential and ground-state wavefunction. The 
length scale for A in / out is given by = h/\/2m±Vo ~ 5 nm, 
assuming a band offset Vo = 20 meV and m± = 0.8mo. With 
these parameters, the asymptotic value of the shell thickness 
(corresponding to the dashed line) is Rout — Rin — 7.7 nm. 


The basic idea is to take advantage of negative band off¬ 
sets and is illustrated in Fig. [l] where we assume now 
that the core (regions 1, 2, 3) and the inner shell (region 
4) each has a distinct chemical composition (e.g., two dif¬ 
ferent SiGe alloys) such that the band offset —Vo of the 
inner shell (region 4) is negative with respect to the core. 
We can estimate a typical thickness of the shell by us¬ 
ing band parameters of electrons in Ge and by assuming 
for simplicity that both the nanowire growth direction 
and band minimum (valley) are along [ 111 ], such that 
the transverse mass is m± — 0.08mo (with mo the free 
electron mass; note that the valley degeneracy is gener¬ 
ally broken by the presence of strain and confinement). 
In this case, the radial motion is described by: 

H ± = - V o 0(r - R in )d(R out - r) + U(r), ( 8 ) 

2m_L 

where i?i n / ou t is the inner/outer radius of region 4 and 
typical values of Vo for SiGe compounds can reach up 
to several hundred meVP^*^ If the outer shell (region 5) 
has a large positive band offset, it can be approximated 
as an infinite barrier 


U(r) 


0 for r < Rout, 
oo for r > .Rout- 


tS) 


The profile of the total radial potential is schemati¬ 
cally represented in the inset of Fig. [2] As in the one¬ 
dimensional problem, the requirement of a uniform wave- 
function in the core corresponds to a zero-energy ground- 
state. Thus, the ground-state wavefunction has the form 
aJo(qr) + /3Yo(qr) in region 4, where J n (x) and Y n (x) are 
Bessel functions of the first and second kind, respectively, 
and q = \j2m±Vo/h 2 . Imposing the relevant boundary 


conditions at both R ln and R ou t gives 

Jo{qRout)Yi(qR in ) - Y 0 (qR out )J 1 (qR in ) = 0. (10) 


For a given core radius R in and chemical composition 
(i.e., the offset Vo), Eq. (10) determines the appropriate 
thickness R out — Ri n - The result is plotted in Fig. [2] as a 
function of qRi n . At qR[ n = 0 the value qR ou t — 2.405 is 
obtained from the first zero of Jq{x). At large values of 
qR ln the asymptotic result is 


QRout QRin — 


( 11 ) 


As a numerical example, Vo = 20 meV gives q~ l ~ 5 
nm (the length scale of Fig. 2]). For this value, Ri n = 
10 nm corresponds to a shell thickness of 8.6 nm. In 
this example, the shell thickness will approach the value 
7.7 nm, from Eq. ©> for larger values of R; n . 

Core/shell structures can have a variety of designs and 
there should be no fundamental limitation to reach a high 
level of accuracy in fabrication (comparable to III-V pla¬ 
nar heterostructures). We will therefore not pursue a 
more specific analysis of an ideal setup, but will sim¬ 
ply assume in the following that a strategy similar to 
that given here should allow for an accurate realization 
of Eq. (|3| . We discuss next how the uniformity of the hy- 
perfine coupling can be revealed through electric trans¬ 
port measurements. 


IV. SPIN VALVE SETUP 

A quantum-dot spin valve consists of a quantum dot in 
contact with two ferromagnetic reservoirs with opposite 
polarization. The isolated dot and leads can be described 
through the Hamiltonian 


#el = £(^~ 


(v-—\dU + l + e (? V r ^ 

\ 9 2 / a ° ~ J y c p(J L pa L pcr i c p a ' pa' P<? J 

pa 


(12) 

where for the quantum dot we have assumed a single or¬ 
bital level is relevant and have included the effect of the 
Zeeman term. The sign of b = —gfj,sB z depends both 
on g and B z but has no effect on our discussion so we 
fix b > 0. In Eq. (12 1 , l pa (r ptr ) destroys an electron 


in the state p of theTeft (right) lead, with a = +(—) 
corresponding to f (I)- The single-particle energies epf^ 
are spin-dependent and determine the densities-of-states 
vi a , v ra at the Fermi levels g,p r . Assuming identical leads, 
we have = v r ± = u + and 1/4 = v r ^ for the ma¬ 

jority and minority carriers, respectively. Furthermore, 
we consider spin-independent tunneling: 

Ht = Y i till a d a + y; t r rl a d a + h.c., (13) 

pa pa 

where h.c. indicates hermitian conjugate terms. 

The full Hamiltonian, including the nuclear bath, reads 


It = //el + Hr + //),f + H 


N, 


(14) 
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FIG. 3. Schematic of the four relevant tunneling pro¬ 
cesses through a quantum-dot spin valve with positive bias 
A fi and large Zeeman splitting b. The polarizations of major¬ 
ity/minority carriers in the ferromagnetic leads are indicated 
with wide/narrow empty arrows. Panels (a) and (d) show 
direct tunneling events, with rates given by Eqs. ( |16| ) and 
(jTTJ). Panels (b) and (c) show second-order tunneling events 
that involve an electron-nuclear spin flip-flop through hyper- 
fine coup ling. The corresponding rates are given in Eqs. (211 
and (22 I. Processes (c) and (d) are suppressed for half-metal 
leads. 


where Hm is given by the two terms in Eq. ^ . The last 
term is an inhomogeneous field acting on each nuclear 
spin: 

Hn = E b k I z ,k, (15) 

i 

which includes the small Zeeman splitting of the nuclear 
spins ( 6 ^ <C 6 ) and accounts phenomenologically for all 
possible sources of inhomogeneity (e.g., spatial variations 
of the magnetic field B, slightly non-uniform couplings 
A k , or the nuclear-spin dipole-dipole interactions). 

In the following, we will be especially interested in 
electron-nuclear flip-flop processes. In fact, transport 
through the spin valve induces a non-trivial nuclear-spin 
dynamics through H^j , defined in Eq. |3|. More specifi¬ 
cally, it is convenient to assume zero temperature for the 
electron system (but not for the nuclear spins) and a rel¬ 
atively large magnetic held, such that only the f state 
of the dot is in the bias window of width A/r = fii — fj, r 
(see Fig. |3j) . For this arrangement, tunneling out of the 
dot is suppressed due to the reduced density-of-states 
for the minority-carrier band (f) in the right lead. In 
the extreme limit of half-metal leads ( v _ = 0 ), direct 
tunneling out of the dot becomes impossible and the 
current to the right lead is dominated by second-order 
processes involving both Ht and Hf f f, as illustrated in 
Fig-! Clearly, each electron passing through the spin 
valve will transfer H of angular momentum to the nuclear- 
spin system. Even with imperfectly-polarized ferromag¬ 
netic leads, such spin-flip processes induce a transfer of 
angular momentum to the nuclear spins. 


We also note that the left ferromagnet in Fig. [3] is 
not strictly necessary for the transfer of spin, when a 
fixed bias direction is chosen. However, an advantage 
of the setup shown in Fig. [3] is that nuclear spins can 
be polarized in either direction by reversing the bias 
(A n —> — Ap,). 


V. TUNNELING RATES AT LARGE 
MAGNETIC FIELD 


We now discuss the general form of the tunneling rates, 
which determine the electron transport through the spin 
valve. It is most transparent to work with sufficiently 
large magnetic field and in the weak tunneling limit, such 
that both Hm and Ht can be treated as perturbations. 
The precise condition on b is given towards the end of 
this Section, together with some considerations on what 
to expect outside the large -6 param eter re gime (based on 
the non-perturbative results of Sec. |VIB ). 

In the large -6 limit, the unperturbed eigenstates are 
IT® 1 ) (g> I’!'™), where IT® 1 ) (|T^)) are the electronic (nu¬ 
clear) eigenstates of H e \ (Hn), with n (m) a suitable 
index. We also restrict ourselves to the case of forward 
bias, A fi > 0, for which the four relevant transport pro¬ 
cesses are illustrated in Fig. [3] The corresponding tun- 

(i ) 1 — (V) 

neling rates are denoted and where we only 

specify the spin orientation (a =t / |) in the right/left 
(r/l) lead, since in the quantum dot ( d ) the spin is al¬ 
ways f- The reverse-bias case, A/i < 0, can be treated 
in a similar way, by considering the relevant rates 
and T« Q . 

The spin-conserving rates due to Ht are (at vanishing 
temperature, T = 0): 

T^ d = 27ru + \ tl \ 2 = T^\ (16) 

rjj£ t = 2nu-\t r \ 2 = rL r) , (17) 

where the weak-tunneling condition requires Tj_’ <C 6 . 
As for the spin-flip tunneling rates, these involve nec¬ 
essarily both Ht and Fjj, thus need to be computed 
with higher-order perturbation theory. By specializing 
to and to an initial nuclear state |TA) (the initial 

electron state is |T|') = dj| 0 )) we have: 


<*£ I <*/" I H t Hf f | T® 1 ) | Tf) 

xS(Ef„+Ey,-Ei), (18) 

where Ei is the eigenvalue of H$ = H e \ + Hn for the 
initial state and (Ej\,) is the eigenvalue of Hn (H e \) 
for the final state. Since Hj j flips the electron spin in 
the dot, we approximate H 0 — E, ~ 6 , which is applicable 
when 6 jv,i <C 6 , thus the change in nuclear energy induced 


r &4. = 27r E 
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by H'jj can be neglected. Then Eq. (18 1 evaluates to: 


r { J\i = EK*Tl / +l*f>| 2 = 

S' 

(19) 

where we have introduced T^ = 2ni' + \t r \ 2 (similarly, 
T { 1 ] = 27rz/_|ti| 2 ) and 


T1 = 


( A 


\2Nb 


( 20 ) 


The other spin-flip rate, Fy/^, can be obtained with 
the same method. The results are immediately extended 
to a nuclear state which is an incoherent mixture of eigen¬ 
states, i.e., with density matrix p^ = JV 
which will be useful in later sections. Finally, we can 
write the spin-flip rates as follows: 


r 


(0 


I^Apn] = V^ { - Tr [I+I-Pn], 

r^[p N } = vr^Tr[i_i + p N }. 


( 21 ) 

( 22 ) 


The validity of Eq. (18) requires the matrix element from 


the initial to the intermediate state to be smaller than 
b. By considering total angular momentum eigenstates 
|/, m) for the nuclear-spin system (where m = —I ,..., / 
is the eigenvalue of J z ), the most restrictive condition is 
obtained for |V/2,0) (assuming spin-1/2 nuclei), giving 
^(Tp |/ + |'F( V ) ~ A/4 <C b. The actual range of mag¬ 
netic fields at which this condition is satisfied depends 
on the specific system, and can be quite accessible for 
group-IV materials. Using A ~ 4.3 peV and g = 2, ap¬ 
propriate for a 29 Si quantum dot,^ gives B z 40 mT. 
This requirement is further relaxed if the number of spin¬ 
carrying nuclei N is only a fraction of N (iot (the total 
number of atoms within a quantum-dot Bohr radius). 
Then, A should be rescaled by a factor N/N^ot < 1. As 
discussed in Sec. |III[ this scenario is relevant to realize a 
uniform hyperfine coupling. 

We also comment about the smallness of rj when b 3> A 
and N 3> 1. To reassure the reader, we anticipate that in 
a suitable regime of long-range nuclear-spin coherence the 
factor Tt[I±I t pn] is of order N 2 . Furthermore, a deriva¬ 
tion of transition rates beyond the perturbative result of 
Eqs. ([21])— (22) is discussed in Sec. VIB and applies when 
the local nuclear fields bk are constant. In this c ase, the 
eigenstates of the quantum dot are known exact IjEmSHH 
and we briefly review this exact solution in Appendix [A] 
We find that the non-perturbative result is qualitatively 
similar to the large-6 limit, but A/b is replaced by a factor 
of order unity [see Eq. (34) and related discussion]. 

A crucial feature of Eqs. (21) (22) is the nuclear-spin 


factor, Tr [I^I±pn\, since this results in a strong depen¬ 
dence of the tunneling rates on coherence properties of 
the nuclear-spin bath. This dependence becomes espe¬ 
cially pronounced when contrasting the regime in which 
the total angular momentum / is conserved with the op¬ 
posite situation, in which a fast local nuclear-spin dephas¬ 
ing mechanism exists. The difference between these two 


limits is illustrated most clearly for the case of half-metal 
leads, which we discuss next. 


VI. THE HALF-METAL LIMIT 

For half-metal leads ( v _ = 0) we assume that a fully- 
polarized nucle ar-sp in state can be prepared. Nuclear- 
spin dark states^ 2 , ®J-i are an obstacle towards reaching full 
polarization but the pure dephasing term Hn can induce 
transitions out of the dark states.^ Thus, in our model, 
this term allows the system to reach the stationary state 
|iV/2,— iV/2) (for spin-1/2 nuclei and A p < 0). After 
initialization, switching the bias to A p > 0 leads to 

+ + 2)(23) 

where we assume nuclear-spin dephasing is sufficiently 
weak that the total angular momentum is conserved dur¬ 
ing the time evolution. In particular, this assumption 
is justified if all the couplings Ak,bk in Eq. (14) are ap¬ 
proximately constant. In Sec. | VI A] below we also address 
the opposite case with strong local nuclear-spin dephas¬ 
ing. The transition rates can be immediately obtained at 
large b from Eqs. (21) (22) but we will also discuss the 
rates for arbitrary b in Sec. VI B[ making use of the exact 
solution given in Appendix |A| 


A. Large magnetic field 

For positive bias, the relevant spin-flip transition rates 
are Eqs. (21) and (22). By making use of the angular- 


momentum states /, m), these rates become: 


= r)T { l\l + m)(I — m + 1), 

m + 1 ). 


.(0 


(24) 

(25) 


This result shows that a large enhancement of the tun¬ 
neling rate can be realized since we assume that full 
polarization can be reached for half-nretal leads. Fur¬ 
ther assuming nuclear-spin /^ = 1/2 for simplicity, we 
have I = N/2 for the collection of N coupled nuclear 
spins. _ The largest enhancement is at m = 0, when 


Eq. (25) gives rj^ d ~ 7 ?r^iV 2 /4 (while j. 
to v _ =0). This enhancement is directly due to quan¬ 
tum coherence in the | N/2, 0) nuclear state. On the other 
hand, the uncorrelated initial state at m = —N/2 gives 
— V^-N. The latter result is proportional to N , 
which is the expected dependence for a tunneling pro¬ 
cess where the electron spin flip-flop with the individual 
nuclear spins occurs incoherently. 

To illustrate this coherent enhancement more clearly, 
we consider an incoherent mixture of ‘product states’ of 
type | n) = \ • • ■) an d fixed value of m. In this case, 

Eqs. (21) and (22) give (for generic imperfectly-polarized 


At) 


= 0 . due 
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ferromagnetic leads with v_ ^ 0), 


B. Rates from exact electron-nuclear spin 
eigenstates 


= t?r^(iv/2 + m), 
r^=r?r«(iV/2-m). 


i(0 


(26) 

(27) 


These rates are appropriate if dephasing mechanisms of 
the nuclear-spin system act quickly on the time scale of 
individual tunneling events. This dephasing will occur, 
e.g., in the presence of strong variations in the b *. of 
Eq. (15). In particular, the relevant rate F^J^ for half¬ 
metal leads is simply proportional to the number of ,1- 
nuclei, which allow the spin-flip tunneling process out of 
the dot. 

The enhancement factor of order N is especially sig¬ 
nificant since a typical quantum dot has N ~ 10 5 — 10 6 . 
However, the effect of nuclear-spin coherence can already 
be seen in the first few individual tunneling events at 


(0 

t—►<* 


d —>4, 5 


the rates for the 


I z ~ —N/2. Since T. 

I z = m —> m + 1 transition of Eq. (23) are approximated 
very well by Eq. (25): 


7m — T]Nr+\ 2 V NT { ;\ 3r)NT+\ 


(28) 


for m + N/2 = 0,1,2,.... These tunneling rates 
should be compared with the approximately constant 
( r ) 

rate t]NT\' for an incoherent nuclear-spin bath [see 
Eq. (27) with to ~ ~N[‘A- Thus, if single tunneling 

> 1 pro¬ 


events can be detectedobserving jm+i/ln 
vides a signature of nuclear-spin coherence. The differ¬ 
ence is largest for the first few tunneling events, e.g., 
(7_«+i)/(7_ iv ) — 1 for an incoherent nuclear-spin mix¬ 
ture and (7_jv + 1 )/(7_n;) ~ 2 in Eq. (28) above (while 
lm+i/lm —> 1 for larger m also in the coherent case). 


The enhancement of spin-flip rates through nuclear- 
spin coherence in such a spin valve setup is analogous to 
the enhancement analyzed in Ref. [52] in the spin-blockade 
regime of a double quantum dot. There, it was noted that 
for uniform hyperfine couplings the transport becomes 
analogous to the sup erra diant emission of an ensemble of 
N two-level atomsPSSU An alternative transport setup 
based on this analogy has also been recently examined 
in Ref. [ST] Thus, by interpreting the f / 4- nuclear-spin 
states as ground/excited atomic states, a spin-flip tun¬ 
neling event corresponds to the emission (or absorption) 
of a photon. Given N initially excited atoms, each with 
an independent decay rate t)T+, Eq. (28) is then im¬ 
mediately understood as the decay rate due to collective 
photon emission.^ The ‘superradiant’ regime can be con- 

(r) 

trasted with the photon emission rate r)T\ ' N± for the 
usual spontaneous emission of IVj. excited atoms (with 
jVj_ = N/2 — m), corresponding to an incoherent nuclear- 
spin bath. 


We can extend the previous results to arbitrary mag¬ 
netic field using the exact eigenstates described in Ap¬ 
pendix |Aj In direct analogy with Eq. (23), and assuming 
that only the lower energy level ej m is in the bias win¬ 
dow, we can write the evolution of the coupled electron- 
nuclear-spin system as 


\I, m)|0) -)■ |<Pj (to)) -+ |/,to+1)|0) ->■- (29) 


The corresponding rates are obtained from Eq. (A2) as: 


n‘^ = Tf\/3 Itr 


TO 
■ j->-d 

rW 
1 d-R 


= r 


(r) 


|^/,m | 


(30) 

(31) 


From Eq. ( |A7| , we immediately recover the results from 
the perturbative treatment at large b, since Eq. (30) 
yields Eq. (16) while Eq. (31) yields Eq. (251. On the 


other hand, these expressions allow us to go beyond the 
perturbative limit and discuss the opposite case, 


Assuming a single energy level, e N / 2 r 
dow, and using Eq. (A5), we obtain: 


TO 


= r 


rM — r 

1 d-k|. — 1 + 


(0 N/2 - to 
N + l ’ 

( r ) N/2 + to -)-1 


N + l 


in the bias win- 

(32) 

(33) 


where, once again, we have assumed the fully-polarized 
initial condition, I = N/2. In this case, the hypothesis 
of a single energy level e)y^ 2 m in the bias window can be 
easily satisfied due to the relatively large energy s plitt ing 


t N/2,r 


— e 


N/2,r, 


~ A/2, which, as seen in Eq. (A6), is 


independent of m when 5 = 0. 

To discuss the enhancement of electronic current 
through the spin valve due to nuclear-spin coherence, we 
compute the rate for the process |1V/2 ,to) —> \N/2,m + 
= (1/r^d + !/ r d r -li) _1 - By assuming 


1), given by 

Tn 

(A (’ r '\ 

for simplicity T 7 = T, , we obtain: 


= r 


( r ) (N/2 — m)(N/2 + to + 1) 
(N + l) 2 ’ 


(34) 


which has the same form as Eq. (25) with the substitu¬ 
tion r) —> (N + l) -2 . Thus, all previous remarks on the 
enhancement of the coherent rates at I z = 0 still hold in 
this case. While decreasing b clearly enhances the pref¬ 
actor r/ = A 2 /(2Nb) 2 in Eq. (251, the perturbative treat¬ 
ment fails at small b and the tunneling rate saturates to 
the value given in Eq. (34). For to = 0, Eq. (341 gives 


~ r^/ 4 , comparable to the direct tunneling rate in 
the absence of any electron-nuclear spin-flip mechanisms. 

























VII. DYNAMICS WITH 
IMPERFECTLY-POLARIZED 
FERROMAGNETIC LEADS 

We now return to the limit of large b and discuss the 
general case of ferromagnetic leads with imperfect polar¬ 
ization. In general, the evolution of the reduced nuclear- 
spin density matrix pn is described by a quantum master 
equation. This approach will be described in Sec. VIII 
where we will also justify in more detail how, in the 
two limiting regimes of very fast/slow nuclear-spin de¬ 
coherence, the nuclear-spin dynamics can be character¬ 
ized through much simpler rate equations. As discussed 
here, the final result can be derived more physically di¬ 
rectly from the transition rates discussed in the previous 
section. 

We start from the incoherent evolution, for which we 
suppose a nuclear-spin state of the form 

pn (t) = (*) =J2 pm (ftp™ ’ ( 35 ) 

m rim \N/2-\-mJ m 


seen above, here p m is maximally mixed in the I z = m 
subspace, which is justified in the presence of fast lo¬ 
cal nuclear-spin dephasing. Transitions m —> m ± 1 in 
the nuclear-spin system are induced by electron tunnel¬ 
ing events with the rates obtained in Sec. |Vj such that 
the p m satisfy: 

Pm = 7m—lPm-l + Im+lPm+l ~ (7m + 7 m)Pm- (36) 

To find the value of yjjj, we use the fact that the spin-flip 
tunneling rates are small: <C r^ d ,T^^. It 

then follows that, at any given moment, the quantum-dot 
occupation n d is approximately determined by the direct 
tunneling processes: 


r 


(0 


n d ^ 




+ r(_ r) ’ 


(37) 


Given the probability n d that the dot is full, a nuclear- 
spin-flip process I z = m —► m + 1 can occur due to a 
flip-flop tunneling event from the quantum dot to the 
right lead: 


7m = n d r d^lPm] = ?7r + (iV/2 - to), 


(38) 


while the rate for to to — 1 is determined by a flip-flop 
tunneling event from the left lead to the empty quantum 
dot: 


7m = (! - n d) Ty 


(0 


,[pm] =vT_(N/2 + m). (39) 


Equations (38) and (39) are computed here using 
Eqs. (21) and (22) with the fully mixed states p m given 


in Eq. (35). We have also introduced: 


(0 r M 


r+ = 


tvt 


r (0 + r M' 


(40) 


We now consider the limit of negligible nuclear-spin 
dephasing, for which a suitable choice of pn is in terms 
of total angular-momentum eigenstates 


PN{t) = ^2pi t m(t)\I,m)(I,m\ 


(41) 


For simplicity, in Eq. (41) we have omitted a sum over the 


permutation indexpSl which does not enter the transition 
rates. The degeneracies will be accounted for appropri¬ 
ately when needed. Because of the uniform hyperfine 
interaction, the flip-flop tunneling processes conserve / 
but induce nuclear spin-flip transitions to —» m± 1. The 


relevant tunneling rates are obtained as in Eqs. (38) and 
(39) but with states \I,m)(I,m\ instead of p m , 


7/m 


= r)T±(I =F rn){I ± m + 1), (coherent) (42) 


The time evolution of pp m follows from Eq. (42) as 


(35) 

Pl,m — 7/,m— lPl,m- 1 + 7/,m+lW,m+l 


As 

- (7 t,m + 7 T,m)Pl,m- 

(43) 


As seen in Eq. (43) above, it is simple to recover the limit 
of half-metal leads discussed in the previous section. In 

(r) 

that case, we have T + = T\ , T_ = 0 and maximal 
angu lar momentum / = N/ 2. In t he h alf-metal limit, 
Eq. (42) recovers the result of Eq. (28), 7jE 2 m = 7™ 

and 7 ~ /2j 


= 0. We also note that the probability 
Pi,m is independent of time due to angular- 


P(/) = E 

momentum conservation, which is straightforward to ver¬ 
ify from Eq. (|43| . In addition, we define p m ( t ) as follows 


Pm(t) = ^2pi,m(t) = P(I)Pm\l{t)- (44) 


Equation (44) can be compared more readily to the in¬ 
coherent result of Eq. (361. In Eq. (44) we have intro¬ 


duced the conditional probability p m |7(t) = Pi,m/P{I)- 
Since P(I ) is time-independent in the coherent regime, 
Pm\i(t) obeys the same equation of motion as pi, m (t) = 
Pm\i{t)P(I), see Eq. Q. 


A. Stationary states 

Before discussing dynamics of the nuclear-spin mag¬ 
netization, here we first analyze the stationary solution 
of Eq. ( [36| , which can be found directly from the anal¬ 
ogy with a system of N 2-level atoms discussed at the 
end of Sec. [V] For independent atoms, the probabilities 
of occupying the t / i states are r±/(r + + T_). These 
populations can be interpreted in terms of a fictitious 
spin temperatures T£ at each nuclear site: 

k n Tl =A k ( In ^ ) =^L. 


r. 


In R 


-l 


(45) 
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where is the energy splitting for nucleus k. A possi¬ 
ble choice for A/, is its time-averaged value (i.e., A*, = 
\bk+nd,A/2N\ if A/N is sufficiently small). However, any 
redefinition of A& can be absorbed into the definition of 
T£, the relevant parameter being A k/T£ or, equivalently, 




(46) 


The stationary nuclear polarization is then given in terms 
of a binomial distribution, 


p eq = 
im 


N 


RN/2~ r 


N/2 + mJ (1 + R) N ’ 


(incoherent) (47) 


which satisfies Eq. (36), as can be easily checked. For 


large N , Eq. (47) is a narrow distribution with 


N fl-R 

( " ,) = t(tT R 


(A To 2 ) = 


NR 


(i + Ry 


(48) 


where (to) = m Pm indicates the average with respect 
to the nuclear-spin distribution and Am = to — (to) . 
The coherent rates given in Eq. (42) can also be in¬ 


terpreted in terms of a single fictitious temperature T* 
since we have iJm/lfm-i = T-/T+ = e ~ A / k sT , i nc } e _ 
pendent of / and to. T£ = T* independent of k in the 
homogeneous limit [A* = A in Eq. (45)], required for 
the coherent regime. However, in the coherent regime, 
the physical picture is quite different since Eq. (43) de¬ 


scribes the thermalization of a (21 + l)-level system with 
constant energy-level spacing A. For a given J, the 
equilibrium probability for occupation of state |J, to) is 

a e -( / - m ) A / fc sT* > which gives 


eq (1 ^ R)R 

Pm\I = 


I-r 


1 - R 2I +' 


(coherent) 


(49) 


Equation (49) is very narrow and almost fully polarized 
within a subspace of fixed /, with (in) — I and (Am 2 ) — 
R /(1 — R) 2 . This distribution is thus very different from 
Eq. (47). However, in general we should also account for 


the distribution of angular momentum P(J), see Eq. (441. 


As a first example, we consider a quantum dot discon¬ 
nected from the ferromagnetic reservoirs, in which the 
nuclear-spin system has relaxed to a fully mixed state. 
Since the number of states with the same value of / and 
to is D(N, I) = ( N % +I ) - (jv/a+j+i), we have: 


P{I) = 2~ N D(N,I)(2I + 1). 


(50) 


By allowing transport through the spin valve (t^ r ^ 0) 
and assuming coherent nuclear-spin evolution, Eqs. (44) 


and (49) give the resulting stationary state. The limiting 


case R —> 1 (normal leads) results in a maximally mixed 
state, p m —► 2~ N ( N / 2 +m) ■ This state is the same as the 
result for the incoherent distribution, Eq. (47), and is 
represented by the thick curve in Fig. [4j Thus, unpolar¬ 
ized leads generally do not modify this initial distribu¬ 
tion. For R —► 0 the maximum polarization, to = /, is 



FIG. 4. Stationary distributions of p m , with N = 500 and 
R = 0.1. The thicker curve (with (m) = 0) is p m for a 
fully-mixed nuclear-spin state, assumed to be the initial state 
of the coherent evolution. The thinner dashed curve is the 
stationary state in the coherent regime. The thinner solid 
curve is the incoherent stationary state given by Eq. \Pl\ and 
is actually independent of the initial state. The inset shows 
the dependence of the stationary polarization p = 2 (to) /N on 
1 — R in the coherent (solid) and incoherent (dashed) regimes. 


reached in each subspace of fixed I. In this case, Eq. (44) 


evaluates to p m >o = P(m) and p m <o = 0. As is made 
clear from Fig. [4j the two stationary distributions (co¬ 
herent/incoherent) with the same R are generally very 
different and, in particular, the nuclear polarization is 
much smaller for the coherent evolution. This difference 
in polarization is because, in the fully-disordered nuclear- 
spin state, the large majority of states have / ~ 0. Thus, 
in the absence of an efficient relaxation mechanism for /, 
a large nuclear polarization cannot be achieved ! 32 * 81 ! 

An expression can be found for P(I) from the station¬ 
ary state under incoherent evolution given in Eq. (47). 


Since Eq. (35) describes a fully-mixed state within each 
I z = to subspace, the conditional probability is simply 
given by P(I\m) = D(N,I)/( N/ ^ +m ), if I > \m\. Thus, 
in this case we obtain 


P(I) = Y, P(I\rn)p% 

m——I 

_ D(N, /)(1 - R^+^R 1 */ 2 - 1 
~ (1-R)(l + R) N ' 


(51) 


This distribution can be realized if a finite bias ±Ap is 
applied across the spin valve for a sufficiently long time, 
exceeding the nuclear-spin coherence time. We note that 
Pm P e ~m does not affect the distribution of I. Thus, it 
is possible to initialize the system at A p < 0 with the in¬ 
coherent distribution, Eq. ( |51| ), and to examine coherent 
dynamics upon reversing the bias. Evaluating Eq. (44) 


with Eqs. (51) and (49) gives the same p ^ previously 


given in Eq. (47). We thus arrive at the interesting re¬ 


sult that the coherent evolution, combined with this P(I) 
in Eq. (51), leads to the same stationary distribution as 


in the incoherent case. 
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This conclusion should not be very surprising. Physi¬ 
cally, if full decoherence has taken place and the nuclear- 
spin system has reached a stationary state, this state will 
also be stationary for shorter time scales (at which it is 
justified to neglect decoherence). The analogy with a fic¬ 
titious thermalization process at temperature T* offers 
an alternative point of view: while the coherent evolu¬ 
tion leads to a fast thermalization process within each 
subspace with fixed J, it does not allow for thermaliza¬ 
tion between sectors with different I. On the other hand, 
the incoherent evolution leads to full thermalization. At 
long times, thermal equilibrium is established globally 
and within each I sector, such that the final state is sta¬ 
tionary for the coherent evolution as well. 


B. Evolution of the nuclear magnetization 


We turn now to the dynamics of the nuclear magne¬ 
tization. For the incoherent case, an exact solution of 
Eq. 
found: 


which is a generalization of Eq. (47) can be 
N\ [f + (t)] Ar t[f_(t)] Ar r 


Pm ^ \N r J (T++r_) N 

with N-f = N/2 ± m and 

f±(i) = T± =FT 0 e- 7?(r + +r - )t , 


(52) 


(53) 


where To determines the initial condition. Apart from 
the obvious remark that To = 0 or sufficiently large t 
yield back the stationary solution of Eq. (471, it is also 


interesting to consider To = T + — T_. This initial con¬ 
dition is the stationary state for negative bias A/i, and 
thus represents a natural starting point for the evolu¬ 
tion of the nuclear-spin ensemble. The resulting nuclear 
magnetization reads 


t) = 


N 1 — R 
~2 1 + R L 


l _ 2 e -? T+( 1 + fl ) t 


(54) 


In the half-metal limit T + = T and T_ = 0, the coher¬ 
ent dynamics is easily understood in terms of the optical 
analogy to superradiance. In fact, Eq. (43) gives 


with 


Pl,m — 2/1 ,m—lPl ,m—l r Yl : mPl,mi 


7 Pm = rjT{I - in)(I + m + 1). 


(55) 


(56) 


This decay rate in the superradiant regime is very well 
known in the quantum optics community (see Ref. [56] for 
a review). At a given value of J, the analytical solution 
for the distribution governing the nuclear magnetization 


i 56 


2rjTIt 

-xy exp 
m) z 


21(1 + to ) _ 2rirH 


I — m 


, (57) 



m 

FIG. 5. Time evolution of p m in the coherent and incoherent 
regimes. Main panel: coherent time evolution of p m with 
N = 500 and R = 0.25. We plot the distributions at times 
t n = nAt, with lyT+Ai = 0.01 and n = 0,1,...,6 (from 
dark er t o lighter color). At t = 0 the distribution is [see 
Eq. ( |47| ] and is close to p ^ at n = 6. The inset shows the 
same result for incoherent dynamics and a larger time step 
pV+At = 0.4. 


for the initial condition \I,I Z = — /), large I /$> 1, and 
t > 1/(TI). Since half-metal leads would allow full polar¬ 
ization of the nuclear-spin bath, in this case it is justified 
to simply set p I}Tn = Pm\(N/2)fii,N/2, giving a complete 
description of the nuclear-magnetization dynamics. 

The general case of imperfectly-polarized ferromag¬ 
netic leads cannot be mapped exactly to this superra¬ 
diant description since Eq. (43) is not of the form of 


Eq. (55). Furthermore, the initial state is generally not 
|/, —I). Instead, it is necessary to consider a mixture 
of different values of I and mn. We take the initial 
condition to be the stationary state at negative bias, 
Afi < 0. In this case, the initial values of p/, m are 
simply obtained from Eq. (47) as P(I\m)p e 3 m , where 

N 


P(I\m ) = D(N,I)/( N , 2 +m ). We then numerically solve 
the simultaneous equations, Eq. (43). An example of the 
resulting time evolution is shown in Fig. [5]where, in con¬ 
trast with the incoherent dynamics (inset), typical fea¬ 
tures of the superradiant behavior are recognized. These 
features include a much faster dynamics (the timescale is 
shorter by a factor ~ 40) and the broad distribution at 
intermediate times (Am becomes of order N ). To make 
this connection explicit through an analytical treatment, 
we start from the approximation: 


7/,m+lW,m+l 7/, m T/,m — 7l, m Pl,m 7 1 ,m-lPl ,m-lj 

(58) 

which is valid when 7 j m+1 pi : m +1 — 7/ m Pi,m has a weak 
dependence m (e.g., if 7 /, m , pi, m are sufficiently broad 
functions of m). By using Eq. (58), we can rewrite the 
time evolution Eq. (43) in the following form: 


Pl,m ^ (7j( m -l - 2/l t m-l)Pl,m-l - ilt,m ~ 7 1,m)Pl,rn- 

(59) 


For I ± m /$> 1, we can approximate the above rates as 
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time to reach zero polarization is to = ln(3.28/)/(27jr/). 
This translates to 


r]T + t 0 ~ 


(1 + R) ln[1.67V(l - R)/( 1 + R)} 
N{ 1 - R) 2 ! 


(coherent) 

(63) 

which is strongly reduced with N, being multiplied by a 
factor In N/N. The corresponding value of ?;T + to for the 
incoherent dynamics obtained from Eq. (541, 


r]T + t 0 ~ 


In 2 

1 + i?’ 


(incoherent) 


(64) 


is independent of N and can thus become much longer 
than Eq. (63). With the parameters used for Fig. [5] 
Eq. (63) gives r]T + t 0 ~ 0.03, in good agreement with 


the numerical evolution, and Eq. ( |64| ) gives ??r + t 0 — 0.6. 
For a more detailed comparison of the two regimes, and 
of the superradiant approximation, see Fig. [ 6 ] 


C. Current dynamics 

The interesting behavior of the magnetization dynam¬ 
ics is accessible through the electron current J. Neglect¬ 
ing the spin-flip contributions to the current, we obtain 
the lowest-order result due to sequential tunneling: 


FIG. 6. Panel (a): Comparison of the numerical results for 
the coherent time evolution of p m (solid) and the superradiant 
approximation, Eq. (57) (dashed). 

0.015, 0.025, 


The values of r]T+t = 


and 0.035 are indicated for each of the three 
curves. Panel (b): plot of the average nuclear magnetization 
as function of rjT+t for the coherent (solid) and incoherent 
(dot-dashed) evolution. The dashed curve is the approximate 
result obtained from the superradiant distribution, Eq. (57). 
For both plots we have used N = 500 and R = 0.25, as in 

Fig. [5| 


7 7 TO — r]T±(I 2 — to 2 ) or, with the same accuracy, 


Jn = 


p(0p( r ) 

~'(0 p( r 


(65) 


Equation (65) gives the leading contribution to the con¬ 


stant background current through the spin valve, inde¬ 
pendent of the nuclear magnetization. On the other 
hand, the correction SJ = J — Jo depends on the magne¬ 
tization dynamics. This distinction is especially clear for 
half-metal leads, when J 0 = 0 and each tunneling event 
through the spin valve is associated with a nuclear-spin 
flip, giving SJ = d(m)/dt. 

For imperfectly-polarized ferromagnetic leads, the to¬ 
tal current J is given by a formula similar to Eq. (651, 


7 r(m — 7 /, m — ? ?(F+ — r_)(7 — m)(J + to + 1 ), (60) but T® and T)!'' are replaced by the total rates 




vr 


which leads to a dynamics of the same form of Eqs. (55) 
and (56 1 and allows us to identify 


i- 


[/jjv] and rL r) +T i j\ i [p N ]. Such an expression for J 


gives the following lowest-order correction to Eq. (65): 


r~r4--r_. 


(61) 


n 2 d v( d\i{p N } 


To obtain the appropriate value of /, we observe that for 
an incoherent initial state p[ l q m , Eq. (51) gives a narrow 
distribution P(I) centered around 


SJ ~ (1 - rid) r+_>. d [pjv] -r 
= (1 -n d ){j~) +n d ( 7 +), 


( 66 ) 

(67) 


where in the second line we used Eqs. (38) and (39) 


(analogous expressions with (y±) -7 (7 f m ) apply to the 


(/) ~ -(to) ~ 


TV 1 - R 
T 1 + R’ 


(62) 


We can thus assume that the initial state approximates 
1 1,—I), with I as in Eq. (62). 


The discussion above shows that Eq. (57) can be gen- 


coherent case), nd is given by Eq. ( |37j , i.e., neglects 
spin-flip corrections to the dot occupation, which is ap¬ 
propriate for this lowest-order expression for SJ (in the 
spin-flip rates). To make a more direct connection to 
the half-metal limit (when SJ = d(m)/ dt ), we can use 
d(m)/dt = ( 7 +) — ( 7 “) and rewrite Eq. (67) as: 


erally used to describe the nuclear-magnetization dynam¬ 
ics, together with Eqs. (61) and (62). In particular, the 


r T dim) , 

5J ~ nd dt + ^' 


(68) 
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coherent case can be approximately described with: 

SJ ~ ( n d + -^ 77 ^- (coherent) (70) 


1 — R } dt 


This formula is obtained from Eq. ( 68 ) using d(m)/dt = 
and ~ i?y+ m (this is justified when the 
rates with I |m| are most relevant and we can approxi¬ 
mate Eq. |42j) witli 7 ^ m ~ ^r ±(/ 2 —to 2 )). As seen in the 
inset of FigTH] Eq. (701 is in excellent agreement with the 
numerical evaluation. A further approximation can be 


realized by using the superradiant distribution Eq. (57) 


FIG. 7. The main panel shows a comparison between the 
coherent (solid) and incoherent (dashed) dynamics of SJ (no¬ 
tice the logarithmic scale of the SJ axis). The two curves have 
the same initial and asymptotic values {SJ —> r]NT-/(l + R) 
at large times). The inset shows a comparison of the coher¬ 
ent result, Eq. (|68|) (dots), to the approximate expression, 


to evaluate d(m)/dt. This gives the dashed curve in the 
inset of Fig. [7] and allows us to find analytic expressions 
for SJ. In particular, the maximum value of the current 
can be estimated as follows: 


SJmFLic — 0.2 


(i ~ R ) 3 

(1 + R) 2 


R 


n d 


1 -R 


r]N 2 T 


+ > 


(71) 


Eq. (70). The solid line is obtained from Eq. (70f using the and Qccurs at the same time g i ven by Eq. pk. Thus, 


le magnetiza- 


exact magnetization, while for the dashed line t 
tion is obtained from Eq. 157) [see also the solid and dashed 
curve in Fig. [6jb)] . In this plot N = 500, R = 0.25, and 
\tl/t r \ = 1 . 


Imperfectly-polarized ferromagnetic leads introduce the 
multiplicative factor n d < 1 (instead of n d = 1). They 
also result in a finite depolarization rate ( 7 “) of the 
nuclear-spin bath, due to spin-flip tunneling from the left 
contact (which is absent for half-metal leads). 

For the incoherent evolution, we can compute an ex- 

( to ) 


plicit result from Eq. ( 68 ) by using Eq. (|54j) for 
Eq. (39) for 7 “. We o: 


and 


stain: 


SJ ~ 


n d - 


R \ d(m) r]NT_ 


with 


1 + R 

d{m) 
dt 


(incoherent) 


dt 1 + R' 

= r]NT + {l — i?)e -I?r+ ( 1+ ' R)4 , (69) 


where the second term (the background contribution to 
SJ) originates from ( 7 “). We also notice that the pref¬ 
actor of d(m)/dt is different from n dl due a transient 
contribution from ( 7 ”). Since d{m)/dt has a simple 
exponential decay, the sign of n d — R/(l + R) deter¬ 
mines if the transient current is a decreasing or increas¬ 
ing function of time. Both scenarios are possible: since 
n d = (1 + Itr/ti^VR)- 1 , the negative sign is realized if 
\t r /ti \ > R~ 3 / 4 > 1 . 

The slow monotonic dependence of the current in the 
incoherent case should be contrasted with the transient 
current peak of the coherent case, shown in Fig. [7] We 
have discussed in Sec. |VII A| that the stationary nu¬ 
clear state pm is the same for the two cases (coher¬ 
ent/incoherent). From Eq. ( 66 ), one concludes that the 


as in the case of the magnetization dynamics, the time 
scale of the coherent current pulse is much shorter than 
for the incoherent dynamics. The current peak has an 
enhancement factor of order N relative to the incoher¬ 
ent case. This factor is clear comparing Eq. (0 with 
the typical incoherent value ~ 77 NT + . These features are 
analogous to the superradiant emission of light, which 
is also in the form of a short and intense pulse. At 
short times, while the incoherent current is essentially 
constant, the coherent evolution shows an exponential 
increase SJ oc e 2r/t , as also typical for superradiant light 
emission . 5 


VIII. INTERMEDIATE NUCLEAR-SPIN 
DEPHASING 

In the previous Sections, we have treated two extreme 
limits: the coherent and incoherent regime in which the 
nuclear-spin dephasing due to the inhomogeneous term 
Hm in Eq. (14) is either absent or very strong, respec¬ 
tively. We wish now to consider a small but finite de¬ 
gree of dephasing, to establish the robustness of the 
superradiant-like behavior: the analysis will thus char¬ 
acterize the relevant timescale for nuclear-spin dephas¬ 
ing, below which the superradiant-like transport can be 
observed. 


asymptotic value of the current is the same as that given 
in Eq. (69). Additionally, at t = 0 the current is the 
same in the two cases (see Fig. [7| but at intermediate 
times the time evolution is dramatically different. The 


A. Crossover timescale 

The question of the relevant time for superradiant-like 
transport is important since the enhancement of tun¬ 
neling rates depends on the coherence properties of the 
nuclear-spin system and the nuclear bath is comprised 
of a large population of spins, up to N ~ 10 5 — 10 6 . 
As is well known, entangled states with a large number 
of particles have short coherence times, which can scale 
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like 1/N (or even 1/N 2 in the case of spatially correlated 
phase noisd^j). However, we find no unfavorable scaling 
with N in this case. 

To understand the final result in simple terms, we re¬ 
fer to the expressions for the spin-flip tunneling rates, 
Eqs. (21) and (22). For an angular-momentum eigenstate 
| N/2, to), all the\n m \pM\n' m ) have the same value but are 
affected very differently by the inhomogeneous broaden¬ 
ing A6: high-order coherences (n m \pN\n' m ) , where \n m ), 
| n' m ) differ on a large number of nuclear spins, typically 
decay with a large rate ~ NAb. However, Eq. (22) de¬ 


pends on pn through the factor Tr[J + J_/jjv], which only 
involves the second-order coherences, i.e., where \n m ), 

| n' m ) differ by a single flip-flop process. This suggests 
that the relevant time scale to observe the superradiant- 
like transport is given by 


t 4> ~ ( A ^) 


-1 


(72) 


This argument indicates that does not scale with N 
and is of the same order as the dephasing time for a single 
nuclear spin (nuclear-spin coherence times of ~ 1 ms have 
been reported in quantum dotsj® and up to ~ 3 hours 
for ionized donors in silicon®), thus that nuclear-spin de¬ 
coherence does not pose a severe limitation to observing 
the ‘superradiant’ transport regime with large N. We 
will confirm this result below with a more sophisticated 
calculation. 


B. Nuclear-spin master equation 

To derive a master equation for the nuclear-spin dy¬ 
namics, we find it simpler to start from a transformed 
Hamiltonian, obtained by applying: 

U = e ^ s+I -- s - I+ \ (73) 

which eliminates the flip-flop terms of the hyperfine cou¬ 
pling Hff to lowest order in A/b. We obtain: 

UHU t ~ (H e] + H T + H N ) + H,, (74) 

where Hi = H zz + SHt is the interaction Hamiltonian. 
5Ht describes the following flip-flop tunneling processes: 


We now consider a natural partition of the Hamilto¬ 


nian, Eq. (74), into electronic and nuclear-spin degrees 


of freedom, with H 0 = H 0 \ + Ht + Hn and the in¬ 
teraction Hamiltonian Hj. Since the electron dynam¬ 
ics are generally much faster than the nuclear-spin dy¬ 
namics, we derive an approximate master equation for 
the nuclear-spin bath starti ng from the standard Born- 
Markov approximation!®® 

Pn( t) = - iTr e i[H7(t),pjv(0) ® Pel] 

/»00 

— Tr e j / dr[H/(t), [H T (t - r),p N (t) <g> p e i]], 

Jo 

(76) 

where Tr e i [...] indicates a trace with respect to the elec¬ 
tronic degrees of freedom and 0(t ) are operators in the 
interaction picture. In particular: 


I±{t) = ^2h,±e 


it ibut 


(77) 


Equation (76) assumes a factorized density matrix p(t) ~ 
Pn( t) <S> Pel) where p e i is the stationary electron state. 
Taking into account the tunneling process (Ht is in¬ 
cluded in H 0 ): 


p e i = n d dj|0)(0|d t + (1 - n d )|0)(0|, 


(78) 


where nd is given by Eq. (|37j). Working in a weak- 
tunneling regime, we have neglected the effect of tun¬ 
neling on the electronic states |0),dt|0) appearing in 


Eq. (78) and, within the same range of validity, we ne¬ 


glect Ht in the interaction picture operators. As a result, 

lpi(t) = Ipie , r„i(t) = r p ±e , which allow one 
to evaluate Eq. (78) in a straightforward waj!®®. We 
also set d-j-(f) = df, by choosing V g = 6/2 in Eq. ( |12[ ). 

A further simplification in the final form is achieved 
with pn{J) diagonal in the I z quantum number. This 
requirement is physically justified since it is consistent 
with the final form of the nuclear-spin master equation 
and, as discussed in Sec. VII A| we can safely assume 
that in the initial state any coherence in I z has decayed 
to zero. Under these assumptions, and returning to the 
Schrodinger picture: 


SHt = ~2M + trV lO d T /+ + hx -’ (75 ) 

P 


arising from the original tunneling Hamiltonian, 
UH t W ~ H t + 5H t . In H u in addition to keeping 
only contributions up to first order in A/6, 6& (we sup¬ 
pose bk ~ A/N ), we have also omitted terms involving 
d® in SHt■ These terms have a negligible influence in 


the regime considered here since df^|0) is outside the 
(large) bias window. Using Eq. ([75]), it is sufficient to 
consider the lowest-order Fermi’s golden rule to obtain 
spin-flip tunneling rates in agreement with Eqs. (21) and 
(| 22 >. 


PN(t) = -i[H N + H L s, PN(t)] + rjr±V[I±\PN{t), 

(79) 

with the standard Lindblad dissipator V[A]p = ApA 7 — 
i(A^Ap + pA 7 A)® and the Lamb-shift Hamiltonian: 

Hls = ^V^lsI 2 , (80) 

where, assuming a uniform density of states within a sym¬ 
metric bandwidth |e®| < A Q : 


p(0 

r LS = — In 

7r 


Pi 

rW 

+ + In 

Hr 

Ai 

7r 

A r 


(81) 
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The final form of the unitary dynamics in Eq. (79) can 
be understood by noticing that, working with a pw diag¬ 
onal with respect to I z , the contribution from H/ { z is ab¬ 
sent. The contribution of SHr{t ) is also zero, due to the 
average over the electronic state. As for the dissipator, 
since Eq. (751 does not conserve the occupation number 
of the quantum dot, all terms in Eq. (76) that are linear in 
SHT{t) vanish. This implies that the two perturbations 
SHt and H zz act independently on the time evolution, 
with no mixed term appearing on the right-hand side of 
Eq. (76 1 . If [ I z ,PN(t )] = 0, the only nonvanishing contri¬ 
bution from Hi(t) is due to SHrit), which yields the final 
result given in Eq. (|79|). As seen, Eq. (79) is consistent 


with our assumption^, pAr(t)] = 0 since both H^ and 
Hls commute with I z and the dissipators T>[I±] preserve 
the diagonal form of pw as well. 

As a final remark on Eq. (79), we discuss the Lamb- 
shift term Eq. (801. In the derivation we have neglected 
small non-uniform terms and used again the fact that 
I±L F can be replaced by I 2 , if [I z , f>N{t)] = 0. A few more 
details are given in Appendix [B] The logarithmic diver¬ 
gence in Eq. (79) at large bandwidth would be cut-off by 
treating the time evolution beyond the Born-Markov ap¬ 
proximation of Eq. (81). However, we refrain ourselves 
from a full microscopic derivation and consider Tig as 
a phenomenological parameter. This is also justified be¬ 
cause terms similar to Hls would appear by including 
higher orders in A/b in the rotated Hamiltonian, Eq. (741. 


C. Numerical results 


We now solve Eq. numerically, for values of b & 
chosen from a Gaussian distribution with standard devi¬ 
ation A b. By focusing on the half-metal regime (T_ = 0), 
an additional simplification arises in the numerical solu¬ 
tion: since the equations for ( n m \pN\n' m ) only depend on 
(n m -i\pN\n' m _i), the problem can be solved iteratively 

for m + N/2 = 0, 1 , 2,_ We assume a fully polarized 

nuclear state at t = 0 , which gives p_ n/2(/) = e~ r]T+Nt . 

Figure [ 8 ] shows an example of the current through 
the quantum dot for different values of A b (due to the 
half-metal leads, the electron current can be obtained as 


J = d(I z )/dt). As expected from Eq. (72), a small/large 
value of A b allows to recover the coherent/incoherent so¬ 
lutions discussed in detail in the previous sections. In 
particular, it is shown that below a timescale t = |A N 1 
(dots), appropriate for the decay of the second-order co¬ 
herence, all curves are well approximated by the coher¬ 
ent evolution. In Figure [ 8 j the deviation of the current 
from the incoherent limit (a simple exponential decay) 
can be easily identified, which demonstrates how the ef¬ 
fect of nuclear-spin coherence is already evident in small 
nuclear-spin clusters (in this example, 8 nuclear spins). 

Of course, the master equation yields much more in¬ 
formation on the statistical properties of the system dy¬ 
namics than the total current. We focus now on the tun¬ 
neling rates to confirm the behavior obtained perturba¬ 



FIG. 8 . Electron current through a quantum dot in contact 
with half-metal leads and 8 nuclear spins. The upper red 
dashed curve is for Ab = 0. Solid curves with Ab/(r;T+) = 
1,2,5,10,20,50 progressively approach the A b —» oo result 
(lower green dashed curve). The dots mark t = /A6 _1 for 
each of the A b values. Tls = 0 and each solid curve is an 
average over 20 realizations of the local nuclear fields. 


tively in Eqs. (21) and (22). If individual tunneling events 
are detected, the observer can keep track of the nuclear 
spin polarization at time t in each single run (with half¬ 
metal leads and starting from | N/2, —N/2)). Under this 
scenario, we introduce the time-dependent tunnel rates 
7 m (f), i.e., the tunnel rates from source to drain, condi¬ 
tional on having nuclear polarization m at time t. The 
7 m (t) depend on the instantaneous nuclear state pN(t) 
and, based on Eqs. (21) and (22), it is natural to expect 


the following behavior: 


m = / ^ n /2 ,m for * « 

l 7m for t > T 0 , 


(82) 


where 7 // 2 m , 7 + are defined in Eqs. (42) and (38), re¬ 
spectively. In other words, at short/long times, the in¬ 
stantaneous tunnel rates yield the coherent/incoherent 
results of earlier sections. We can obtain the (t) by 
solving the detailed balance relation: 

Pm(t) — 7m—1 (t)Pm—1 (t) 7m(t)Pm(t), (83) 


where p m (t) are obtained from Eq. (79). As seen in Fig. |9j 
the 7 m (t) satisfy Eq. (82). At intermediate times, partial 
dephasing of pn entering Eqs. (21) and (22) interpolates 
between the two results. The numerical evaluation of 
7 m (t) in Fig.[9]shows quite clearly the crossover between 
the two regimes, and that all the 7 m (f) approach the 
incoherent result around t ~ 77 , ~ A 5 _1 . 

While in Figs. [ 8 ] and [9] we assumed for simplicity 
Tz,s = 0, we consider now the effect of a finite Lamb-shift. 
It was already discussed that terms oc I 2 , see Eq. (80), 
make the coherent evolution more robust t 5 7 This bel 


nav- 


ior is demonstrated in Fig. (10) where larger values of 


T^s modify the current dynamics and enhance the cur¬ 
rent peak, which approaches the coherent superradiant¬ 
like result. 



























15 



FIG. 9. Instantaneous transition rates 7 m {t) for a system 
with half-metal leads and IV = 10 nuclear spins (initial nuclear 
state is 15, —5}) and Tls = 0. Solid lines are for riT^/Ab = 0.5 
and show a relatively small decay from the t = 0 coherent 
rates ??r+(5 — m )(6 + m). The dashed lines are for rjT+/ Ab = 
5 x 10 -3 . On the 77 , ~ A 6 _1 timescale, they show full decay 
to the incoherent results ryF+(5 — m). 



FIG. 10. Effect of the Lamb shift, Eq. (801, for a sys¬ 
tem of IV = 8 nuclear spins in contact with lialf-metal leads 
and Ab/(r]T+) = 5. The solid curves are with Fns/(r;r+) = 
0, 2,4, 5, approaching the A& = 0 result (upper dashed curve) 
for larger values of Vls/{i 7r +). The lower dashed curve is the 
incoherent result (A 6 = 00) 


the polarization of the leads. By inverting the bias of 
the spin valve, the sign of nuclear polarization can be in¬ 
verted through a fast coherent dynamics. The transport 
current in this coherent regime reveals features analogous 
to the superradiant light emission, ^ as co nsidered the¬ 
oretically in several recent works. 51 52 57 l 6 ° l In particular, 
the enhancement of the transient current by a large fac¬ 
tor of order ~ N (the number of nuclear spins) is due to 
the creation of long-range coherence in the nuclear-spin 
system, which thus can be directly monitored through 
electron transport. 

We also analyzed the crossover between the two (coher¬ 
ent/incolrernt) regimes, showing that the collective en¬ 
hancement survives on a relatively long timescale, i.e., 
the dephasing time of individual nuclear spins. Further¬ 
more, it is well known that the largest collective enhance¬ 
ment occurs in the limit of uniform coupling constants. 
We have outlined a realistic strategy to realize this limit 
in nanowire quantum dots, through the fabrication of a 
“nuclear-spin island” of spinful isotopes, complemented 
by wavefunction engineering in the longitudinal and ra¬ 
dial directions. The latter can be realized through a suit¬ 
ably tailored core-shell structure. We expect that the 
implementation of such “box model” hyperfine interac¬ 
tion, by allowing a controlled interaction of the electron 
spin to the nuclear-spin system, could find many interest¬ 
ing applications well beyond the transport model studied 
here. 
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Appendix A: Exact hyperfine eigenstates with 
uniform coupling 


IX. CONCLUSION 

We have discussed transport through a quantum dot 
in contact with ferromagnetic leads, for a spin valve con¬ 
figuration realized by several recent experimentsP^HMI 67 ! 
Focusing on the hyperfine-mediated flip-flop processes, 
we have analyzed two distinct regimes of coher¬ 
ent/incoherent evolution. In the coherent limit the 
nuclear-spin system is quickly driven into dark states 81 
and only a small amount of polarization can be gener¬ 
ated, starting from an unpolarized thermal state. How¬ 
ever, nuclear-spin dephasing can drive the nuclear-spin 
system away from such dark states, resulting in a sizable 
nuclear polarization in the long-time limit. This station¬ 
ary value of the nuclear polarization is simply given by 


We summarize here the exact solution of the quan¬ 
tum dot Hamiltonian with uniform hyperfine coupling, 
7?hf + H c \, see Eqs. ^ and ©. If |0) describes the 
electronic state with an empty dot, the eigenstates are 
simply given by |0)|/,m) (suppressing an additional per¬ 
mutation index, which does not enter the eigenvalues®!), 
where I(I + 1) and to are the eigenvalues of I 2 and I z , 
respectively. For a singly-occupied dot we have the fol¬ 
lowing eigenstates (with to = — I... I — 1): 

\(pj(m)) = a/ im d||0)|/,TO) +/9/, m d||0)|/,TO+ 1), 

(Al) 

\<pj{m)) = a/ jm d[|0)|/,TO + 1) -/3j, m <4|0)|I,TO), 

(A2) 
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where a/ iTTl = cos(0j, m /2) and /3/ im = sin(0/ im /2), with 
6 A 


9i, m = arg 


4N 


(2m + 1) 


A 


+*— v / i(7Tiy Tr m^n+T) 


■ (A3) 


The corresponding energies are: 
A 


lb 2 bA , 

[ A . A 

/ —- —(‘2m + 1) + 

/ 4 4 N 

[4JV (2/ + 1) J 


*I,m = V 9~4x ± 

(A4) 

The states in Eqs. (All and (A2) are supplemented 
by the fully-polarized eigenstates dj|0)|/, I), dj|0)|/, —I), 
with eigenvalues V g =F b/2 + AI/ (2TV). 

An interesting limiting result is when the external mag¬ 
netic field is zero. By setting 6 = 0 and assuming for 
definiteness A > 0, Eqs. (A31 and (A4) yield 


(the applied bias is A/j = p L — p r , see Fig. |3j) and 
'D^kk' are defined as follows (all operators at time t): 

'D > Ck'[PN{t)] = Ik',- f- Pn Ik - - Ik-Ik',+PN > (B2) 

Hk,k'\PN(t)\ = Ik,- Pn Ik’,+ — Pn Ik’,+h-- (B3) 

The terms of Eq. ( |B1[ ) proportional to n a p can be eval¬ 
uated as follows (the remaining terms with 1 — n a , p can 
be treated in the same way): 


E 


=^c4 


/ o 

f 

/—a q 


n a , p e ±l{ ^- bk ’ )T dT 


de (a) 


pi n “,p 


'x&kpi ~ b k ’) ±ip — 0 


p-t 


- bk< 


(B4) 


(A5) 

(A6) 



6 A , 


term in the square parenthesis of Eq. (B4) evaluates to 

± 

-(m-f 1/2) 

2 2 N y ' _ 

, at 6 A, 

±iv a i In 

p. a -by 

A a +b k , 

and contributes to the Lamb shift Hls 


// + m + 1 // — to 

21 + 1 ’ /3^ ’ m = V 27 TT , 

4,m = V 9 ~^ ± ^(I + 1 / 2 X at 6 = 0, 

In the opposite regime of large magnetic field: 

ui,m ~ 2^6^^ + !) - m ( TO + 1)) Pi,m — 1, (A7) 


~ v _ 

/,rn g 


where in the above expressions we have assumed 6 > 0 
and neglected terms of higher order in A/6. 


Appendix B: Lamb shift Hamiltonian 

By explicitly writing the r dependence of the inte¬ 
grand, the second line of Eq. gives: 

^^ poo 

E W dT\t a \ 2 {n d (l-n a , p )viy,[p N {t)} 

a,p,k,k f ^ 


where we have written ]T) p ~ v a i deff, by assum¬ 
ing a symmetric band with constant density of states 
(z/q = and z/ r | = v + ). Since 6*, -C |/x a |, the integra¬ 
tion of £(e/,“ — bk>) in the second line is immediate. It 
gives 7T v_8 a i, which contributes to the Lindblad dissipa- 
tor of Eq. (79|. 


On the other hand, the integration of the second 


appearing in the first term of Eq. (791. The presence of 


b k > makes the result inhomogeneous with respect to the 
nuclear index k'. However, notice that the effect of b k ' on 

~ In \p a /A a \ -b k >/n a and 


Hls is small since In | + b b ' 


bk'/Pa 1. Therefore, the corrections are smaller than 
II HLsbk'/pa\\ < ||//ls||, \\H n \\, and it is justified to ne¬ 
glect them. It then becomes straightforward to evaluate 
H L s from Eq (|B1|: 


Hls =2 ? 7 r iS 


rL r) /+/- 


■T ( 1 ) I_I+ 


~'( r ) 


''(O 


(B5) 


+(1 — nd)n a , p 'D k fc/[piv(/)]|e*^ e p/ bfc, - )T + h.c., (Bl) where Tls is defined in Eq. (81). As discussed in the 


where a € {/, r} labels the two leads. In Eq. (Bl) we have 


introduced the lead occupation numbers n ajP = 9(p, a — in Eq. (801. 


main text, our assu mpt ion on p jy allows us to substitute 
I±Izp -A 1 2 in Eq. (B5), which leads to the result cited 
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